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fundamental concern when 
A possi or specifying the 

performance parameters of, 
rotating machinery with fluid bear- 
ings is, “Will my machine be stable?” or 
alternatively,“ When will my machine 
be UNstable?” The purpose of this 
article is to provide some information 
about stability, specifically, what stabil- 
ity is, how to measure it, and what fac- 
tors contribute to stability or the lack 
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thereof. In the last issue of Orbit, we 
presented a discussion of the synchro- 
nous amplification factor, Os. This arti- 
cle is an extension of that previous 
one, since amplification factor and sta- 
bility are frequently mentioned in the 
same breath, It is important to note 
that the amplification factor is NOT a 
measure of stability. 

Stability is defined as the resistance 
of a system to disturbance. For the 
purposes of this article, stability is 
concerned with the rotor’s resistance 
to lateral motion, given the presence 
of some disturbing force on the rotor. 

Stability can also be defined in 
terms of its negative, instability. The 
condition of instability dictates that 
the amplitude of the rotor's response 
motion will increase without bound 
even without any external disturbing 
force. In reality, it is impossible for 
motion to increase without bound. 
Typically, the rotor motion will 
increase until system nonlinearities 
(for instance, journal contact with the 
bearing surface) restrain the rotor in 
some sort of limit cycle of self-excited 
vibration. 


But within the linear range of the 
system, the forced response motion of 
the rotor model can be summarized 
as: 


ROTOR RESPONSE MOTION = 
FORCE 
STIFFNESS 


where, for a simply supported rotor 
with one lateral mode, the Dynamic 





Stiffness is defined as follows [1]: 
Total Dynamic Stiffness = 
Direct Dynamic Stiffness + 


Quadrature Dynamic Stiffness 


K _Mw* = Direct Dynamic Stiffness 
(DDS) (1) 
jD(@-ALM) = Quadrature Dynamic 
Stiffness ( ODS ) (2) 


There is a resonance when either 
of these terms become zero. When 
the DDS becomes zero, it is mechani- 
cal resonance by definition. The only 
term limiting the amplitude of motion 
at mechanical resonance is the 
Quadrature Dynamic Stiffness. This 
observed resonance peak has general- 
ly been called the “critical.” When ODS 
hecomes zero, it is fluid induced reso- 
nance, The only term limiting the 
peak of the motion at fluid induced 
resonance is the Direct Stiffness term. 
For the typical case of a lightly 
damped machine, the mechanical res- 
onance will be the dominant factor 
determining system response. 

In a simple spring-mass-damper sys- 
tem, the only Quadrature Stiffness 
term is D@ (damping times excitation 
frequency). This is why many books, 
past and present, use the term “damp- 
ing becomes negative” or “damping 
becomes zero.” In fact, it is impossible 
for damping (D) to be negative. It is 
the Quadrature Dynamic Stiffness that 
becomes negative [2]. Once more: 
When you see a book or paper that 
states that the damping goes negative, 
beware. 


June 1996 


Both resonances are always pre- 
sent, regardless of machine parame- 
ters, however, the  fluid-induced 
resonance has been rarely observed, 
because it requires nonsynchronous 
excitation. To illustrate these reso- 
nances isolated from each other, 
Figure l(a) shows a Bode and polar 
plot of system response to nonsyn- 
chronous perturbation (exciting force 
proportional to @F,as with unbalance) 
with QDS artificially held to a con- 
stant, nonzero value (ODS must be 
nonzero or the amplitude at reso- 
nance is infinite). Similarly, Figure 1(b) 
shows a Bode and polar plot of sys- 
tem response to nonsynchronous per- 
turbation with DOS artificially held to 
a constant nonzero value. The com- 
bined system response, including both 
DOS and GDS influence (not con- 
stant), is shown in Figure l(c). The 
machine parameters are the same for 
all cases (Fig, la, b, and c). Please note 
that these are artificial conditions cre- 
ated to illustrate basic behavior. 

From Eq.¢1), the mechanical reso- 
nance frequency Oesppsp OCCUrs 
when DDS = 0, thus 


K 3) 


= Oeo = IET 

Similarly, from Eq. (2), the fluid- 
induced resonance occurs at an exci- 
tation speed that causes ODS to 
become equal to zero: 


O = esros) = Add (4) 


An cigenvalue analysis of the rotor 
model [1] will show that the rotor 
with fluid bearings is unstable if the 
mechanical resonance occurs at a 
lower speed than the fluid-induced 
resonance. This information, com- 
bined with Eqs. (3) and (4), allows us 
to derive an expression for the Margin 
of Stability. The Margin of Stability is 
defined as the perturbation speed dif- 
ference between zeros (roots) of 
Direct and Quadrature Dynamic 
Stiffness as functions of frequency. 
The margin of stability has units of 
frequency (rad/sec). 


June 1996 


(a) 


(b) 


(c) 








f 
p 
= 
g 
: 
i 
z 
Mechanical resonance; shaft speed = 100 rad/sec. 
K -- Mw 2 +7 1872.5 
F 
E 
= 
= 
a 
= 
Z 





im l 
radise 
Fluid induced resonance; shaft speed = 100 rad/‘sec. 


297.75 +jD(w@—- AQ) 


-J15 
a - 
¥ 45 ae 
T s. 
ee a 
E „iË 
relee ' 





AMPLITUDE 
Š = h W b E 


ñ šo ii iSi è 3 
Fadiset 
System resonance: shaft speed = 100 rad,‘sec, 


K -- M@* + jD (@-- AQ) 1872.5 


Parameters 
K=let Me=1 D=35 A®.45 Q= 100 rad/sec 


Figure 1 


Bode and polar plots of isotropic rotor one lateral mode response to 
nonsynchronous excitation 
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CMS =O rese pos) T Ë restons) 
= \* 10 (5) 
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pretation of the Margin of Stability 
with emphasis on the effect of chang- 
ing shaft speed, £2 As £2 increases, the 
zero crossing of QDS moves closer to 
the zero crossing of DDS, thus reduc- 
ing the Margin of Stability. 

In order to be able to compare var- 
ious systems, it is reasonable to intro- 
duce the Nondimensional Margin of 
Stability (WLS) as follows: 


= 
CNMS) = -L a Ne) 
| Dyes DDS) K 


\M 

At the numerator of Eq. (G) there is 
the Margin of Stability; the denomina- 
tor contains the direct resonance fre- 
quency. For stable systems, the value 
may vary between one and zero. A 
negative value means that the system 
is unstable. The Nondimensional 
Margin of Stability better describes 
and quantifies the system stability sit- 
uation than the Margin of Stability. 


Conclusion 

We have discussed the Margin of 
Stability and Nondimensional Margin 
of Stability as a means of evaluating 
stability. It is assumed that the 
machine parameters (A, M, D, A) are 
known. In this case, practical applica- 
tion of the formulas for stability is as 
casy as plugging these values into the 
aforementioned equations. There 
need be no guesswork when evaluat- 
ing stability of rotor systems using 
these methods, 


These parameters are not common- 
ly known, however, they are readily 
obtainable. In fact, if data has been 
collected and stored (as from balanc- 
ing) it may be possible to find these 
parameters from existing data without 
disturbing the machine, If purchasing 
a new machine, it is possible for the 
vendor to supply data from which 
parameters can be calculated. 
References [3] and [4] describe in 
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DDS ( K - mo” ) For Q; and Q2 > Q} 





Direct and quadrature 
dynamic stiffness 


QDS (jD (@-AQ)) 


Plot of Direct and Quadrature Dynamic Stiffness illustrating change in Margin of Stability 
(MS) as rotor speed ( £2) increases. Dashed line corresponds to Ny, with £2) less than £25. 
Note that DDS dos not change. DDS is independent of rotor speed, The zeros of ODS are 


o= AN; and w= Af) respectively. 


The zero of DDS occurs at: 


w=. K /M 
Figure 2 
detail how to go about finding the References: 
machine parameters, and are available L. Muszynska,A.,"One Lateral Mode 


by contacting the Orbit magazine. 

The pioneering work of defining 
the rotor modeling techniques that 
incorporate these parameters, along 
with the techniques for experimental- 
ly determining them, has been pro- 
gressing for many years. The results 
from the work of engineers, such as 
J.M. Stone and A.E Underwood at 
General Motors in 1946, E.H. Hull at 
General Electric in 1959, and, more 
recently, by Donald Bently and Bently 
Rotor Dynamics Research Corpor- 
ation, have led to accurate modeling 
techniques that have proven invalu- 
able tools for machinery diagnostics 
and design. 

This article is part of a series that 
presents discussions of the fundamen- 
tals of quantifving and discussing 
rotor response, in the interest of pro- 
molting more consistent and correct 
communication between technical 
professionals concerned with rotat- 
ing machinery. In the next issue of the 
Orbit, we extend the results of this 
and the previous article, Synchronous 
Amplification Factor Qs [5], by pre- 
senting the more general case of non- 
synchronous excitation, W 
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